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Abstract 

We derive area limit laws for the various symmetry classes of staircase polygons on the 
square lattice, in a uniform ensemble where, for fixed perimeter, each polygon occurs with the 
same probability. This complements a previous study by Leroux and Rassart, where explicit 
expressions for the area and perimeter generating functions of these classes have been derived. 

MSC numbers: 82B41, 05A16, 39A13 

1 Introduction 

Models of planar lattice polygons have a rich history. The most challenging member is self-avoiding 
polygons 30J, a model of interest not only in enumerative combinatorics, but also in the natural 
sciences such as physics and chemistry. Some solvable subclasses, obtained by imposing convexity 
or directedness constraints, have been understood in detail, see [5] for an overview. An important 
example on the square lattice are staircase polygons, also called parallelogram polygons. These have 
been enumerated with respect to various parameters such as perimeter, area and generalisations 
thereof, site-perimeter, and radius of gyration, see [TUl El [3^1 12H] and references therein. For some 
of these parameters, their distribution has been asymptotically analysed. An example is the area 
distribution in a uniform ensemble where, for fixed perimeter, all polygons occur with the same 
probability. We will call this ensemble the uniform fixed perimeter ensemble. For a large class of 
planar polygon models including staircase polygons, the Airy distribution emerges as the limit law 
of area [HI [3Hj ■ There is compelling numerical evidence that the Airy distribution also appears 
as the limit distribution of area in self-avoiding polygons [351 |30] . 

When counting polygons, translated copies of a polygon are identified, but rotated or reflected 
versions of a polygon are often distinguished. Motivated by applications such as benzenoid count- 
ing in chemistry [33], one may be led to identify objects that are related by point symmetries. 
The task arises to enumerate classes of polygons up to symmetries of the underlying lattice. Two 
different counting problems occur. One may ask for the number of polygons which are fixed by a 
given subgroup of the lattice symmetries. One may also count the number of polygon orbits with 
respect to a such a subgroup. These two counting problems are related. In fact, it is possible 
to express the generating function for orbit counts in terms of those for fixed point counts, via 
the Lemma of Burnside. This has been studied for various classes of column-convex polygons by 
Leroux and co-workers [271 1261 121) , and explicit expressions for generating functions have been 
obtained. 

Symmetry classes of staircase polygons on the square lattice are analysed in [26j . by two 
different approaches. One of them uses a bijection between staircase polygons and Dyck paths 
due to Delest and Viennot [9]. By this bijection, some of the polygon symmetry classes can be 
identified with corresponding symmetry classes of Dyck paths. The toolbox of Dyck path analysis 
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can then be used to solve the corresponding polygon counting problems. A second approach uses 
the Temperley-Bousquet-Melou method [5 to obtain explicit expressions for the perimeter and 
area generating functions. In addition, it is shown in [26] that the number of polygons fixed by some 
non-trivial symmetry is asymptotically negligible to the total number of polygons. This implies 
that the number of symmetry orbits is asymptotically equal to the total number of polygons. 

In this article, we derive the area limit laws for the symmetry subclasses of polygons, within 
the uniform fixed perimeter ensemble. Wheras the result for the full class of staircase polygons 
is known already [12\ 138] . the corresponding problem for the various symmetry subclasses has 
apparently not been studied before. We will use an elementary self-contained approach, which 
is based on a simple decomposition of staircase polygons 39 . This yields, for every subclass, a 
g-difference equation for its perimeter and area generating function. Explicit expressions could be 
obtained from this equation, which might differ from those of |26j . We do not focus on solving, 
but on manipulating the g-difference equation in order to derive the limit law, by an application 
of the moment method. Whereas this approach has been used previously in different contexts, 
see e.g. [42j |43l [TT1 [121 E2] , we would like to stress that our alternative derivation, based on the 
method of dominant balance [38, 40 , simplifies the analysis. As a consequence, corrections to the 
asymptotic behaviour might be mechanically obtained, compare |37j . 

Our results extend those of |26j to a complete symmetry analysis, and to a detailed discussion 
of asymptotics. We remark that the area limit laws cannot easily be obtained from the functional 
equations in 26J, or from the explicit expressions given there. One can use the asymptotic de- 
scription of Dyck paths by Brownian excursions [TJ [2] to obtain part of our results. Whereas this 
requires some results from stochastic processes, our approach only needs basic probability theory 
and singularity analysis of generating functions, as e.g. described in [TB] . 

The plan of this paper is as follows. We will discuss the various square lattice symmetry 
subgroups, and characterise the associated polygon classes by certain decompositions. This induces 
functional equations for their generating functions, from which the area limit laws are obtained, by 
an application of the moment method and the method of dominant balance. We will also indicate 
how some of our results may be obtained within a stochastic approach, and we finally discuss 
extensions of our results. 

2 Symmetry classes and functional equations 

We explain the models, introduce basic constructions, and set the notation, following [39] . We will 
then derive functional equations for the perimeter and area generating function of the symmetry 
subclasses of staircase polygons. 

Consider two fully directed walks on the edges of the square lattice (i.e., walks stepping only 
up or right), which both start at the origin and end in the same vertex, but have no other vertex 
and no edge in common. The edge set of such a configuration is called a staircase polygon, if it is 
nonempty. For a given staircase polygon, consider the construction of moving the upper directed 
walk one unit down and one unit to the right. For each walk, remove its first and its last edge. The 
resulting object is a sequence of (horizontal and vertical) edges and staircase polygons, see Figure 
[TJ The unit square yields the empty sequence. It is easy to see that this construction describes a 
combinatorial bijection between the set V of staircase polygons and the set Q of ordered sequences 
of edges and staircase polygons. Let us denote the corresponding map by / : V — > Q. Thus, for 
a staircase polygon P e V, we have f(P) = {Qi, Q n ) £ Q, where Q, is, for i — 1, . . . , n, 
either a single edge or a staircase polygon. We denote the single horizontal edge by eh, and the 
single vertical edge by e v . The image of the unit square is the empty sequence n = 0, which we 
occasionally identify with a single point, denoted by pt. A variant of this construction will be used 
below, in order to derive a functional equation for the generating function of the staircase polygon 
symmetry classes. 

The perimeter of a staircase polygon P 6 V is defined to be the number of its edges. The 
half-perimeter equals the number of (negative) diagonals no(P) plus one. The area n\(P) of a 
staircase polygon P is defined to be the number of its enclosed squares. It equals the sum of the 
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Figure 1: The set of staircase polygons is in one-to-one correspondence with the set of ordered sequences 
of edges and staircase polygons. A corresponding combinatorial bijection is characterised by shifting the 
upper walk of a staircase polygon one unit down and one unit to the right, and by then removing the first 
and the last edge of each walk. 

/ I \ 

(d u d 2 ) (r) (h,v) 

/ I \ I / I \ 
(d x ) (d 2 ) (r 2 ) (h) (v) 

Figure 2: The lattice of subgroups of D4. The rotation about n/2 is denoted by r, the reflections in the 
positive and the negative diagonal are denoted by di and di, and the reflections in the horizontal and 
vertical axes are denoted by h and v. The identity is omitted. 



lengths of its (negative) diagonals. See Figure Q] for an illustration. The weight of a staircase 
polygon P is the monomial wp(x,q) = x n °^ p > +1 q ni ( p > . The half-perimeter and area generating 
function of a subclass C C V of staircase polygons is the (formal) power series 

C(x,q) := ^2 w P (x,q). 
Pec 

Observe that for eh, e v , pt, and for P G V we have [35] 

Wf-i (pt )(x,q) = x 2 q, w f -i {eh) (x,q) = x 3 q 2 , w f -i^(x, q) = x 3 q 2 , 
w f -i { p)(x,q) = x 2 q ■ w P (xq,q). 

For a polygon P 6 V, consider f(P) = {Qi, ■ ■ ■ , Q n )- In order to retrieve P from (Qi, . . . , Q n ), 
translate Pi = f~ 1 (Qi) in such a way that its lower left square coincides with the upper right 
square of Pi-i, for i G {2, . . . , n}. We say that P is the concatenation of (Pi, . . . , P n ), and write 
P = c(Pi, . . . , P n ). The weight wp(x, q) of P is retrieved from the weights of Pi, . . . , P ra via 

w c (p u ...,p n )(x,q) = 1 w Pl (x,g) • . . . • wp n (x,q). 

Denote by V C V the subset of polygons P = / _1 (P), where P € PU{ei, e u }. We have established 
a combinatorial bijection between the set V and the set of ordered sequences from P. 

The group of point symmetries of the square lattice is the dihedral group ©4. Its non-trivial 
subgroups are depicted in Figured) Note that the above decomposition respects any subgroup of 
the square lattice point symmetries. This observation is the key to deriving functional equations 
for the generating functions of the symmetry subclasses. In the proof of the following proposition, 
will treat two cases in some detail, the remaining ones being handled similarly. 

Proposition 2.1. The half-perimeter and area generating functions of the staircase polygon sym- 
metry subclasses satisfy the following functional equations. 
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1. Class V of all staircase polygons with generating function P(x 1 q): 

2 

P(x, q) = X -^—, r. (2.1) 

v ,q) 1 - 2xq - P(xq, q) v ' 

2. Class S of (r 2 ) -symmetric staircase polygons with generating function S(x, q): 

S(x, q) = 4- (1 + 2xq + S(xq, q)) P(x 2 ,q 2 ). (2.2) 
x A q 

3. Class of (di) -symmetric staircase polygons with generating function D\(x,q): 

Di{x,q) = — r. 

1 - D^xq^q) 

4- Class of (d 2 ) -symmetric staircase polygons with generating function D 2 (x,q): 

D 2 (x, q) = 4- (1 + D 2 (xq, q)) P{x 2 ,q 2 ). 
x z q 

5. Class of (di,d 2 )- symmetric staircase polygons with generating function D\ 2 (x,q): 

1 - — - ^ — > 2 „2\ 



D h2 (x, q) = — (l + D 1>2 (xq, q)) D x {x\ q 2 ). 
x z q 

6. Classes of {h)~, (v)-, and (h, v)-symmetric staircase polygons with generating function H(x, q): 

2 2 ^ " r n 

H(x, q) — x qH{xq, q) + x q- 



1 — xq 

7. Classes (r) -symmetric staircase polygons with generating function R(x,q): 

R(x, q) = x 2 qR{xq, q) + x 2 q. 

Proof. Denote the induced group action a : V4 x V — > V by a(g, P) = gP. 

1. The bijection described above implies the following chain of equalities, compare also 



P( x ,<l) = Yl E w c (Pi,...,F n )(a:,?) 



n =°(p 1 ,...,p n )£(vy 



= E^ E 



w Pl (x,q) w Pn (x,q) 



x 2 q x 2 q 
™=° (p 1 ,...,p n )e(P) n 



x 2 q 



E (^E^,9)) 



(2.3) 



x 2 q- 



lk ( W f-Heu) ?) + W f-He„)(x, q) + E pe p w f -i {P) (x, qfj 



9 

x q 



1 — 2xg — P{xq, q) 



2. For P G 5, we have /(P) = (P x , . . . , P„, C, r 2 P„, . . . , r 2 Pi), where C G S U {e v ,e h ,pt}, an d 
Pi G V U {e v ,eh} for i = 1, . . . , n, compare Figure [H In analogy to the definition of V above, 
define S G V as the pre-image of 5 U {e„, eh,pt} under /. Note that concatenation of Q G V with 
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Figure 3: r 2 -symmetric polygon and corresponding sequence of polygons and edges 



the unit square results in Q again, and that we have WQ(x,q) k — WQ(x k ,q k ). With P(x,q) as 
above, this yields 



S(x,q) = 



E 



w c(P 1 ,...,P n ,C,r 2 P rl ,...,r 2 P 1 )( x , <l) 



1 (Fi,...,P„,C)e(7')"xS 



ces 



n=0 



(F 1 ,...,P n )G(P)™ 



wpi{x,qY 

(x 2 q) 2 



(x 2 q) 2 



(x 2 q + 2x 3 q 2 + x 2 qS{xq, q)) 1 

1 — [2x 2 q z + P [x z q z , q z )) 

1 + 2xq + S(xq, q) 



(2.4) 



4 2 



x 2 q 



l-2x 2 q 2 - P(x 2 q 2 ,q 2 ) 



= —(l + 2xq + S(xq,q))P(x 2 ,q 2 ) 
x £ q 



where the sum over Q € P in the third equation is treated as in eqn. (12.31) . In the last step, we 
applied eqn. (|2.ip . 

3. For a (<ii) -symmetric polygon Q, we have /(Q) = (Pi, . . .P n ), with a (di)-symmetric Pi for 
i = 1, . . . ,n. A calculation similar to that in eqn. (|2.3p then yields the assertion. 

4. For a (c^-symmetric polygon Q, we have f(Q) = (Pi, . . . P„, C, c? 2 P n , . . . , d 2 Pi), where Pi G 
P U {e/t, e„} for i = 1, . . . , n, and where C is a pt or (e^-symmetric. Now a computation similar 
to that in eqn. (|2.4p yields the assertion. 

5. For a (di, d%) -symmetric polygon Q, we have /(Q) = (Pi, . . . P n , C, d 2 P n , ■ ■ ■ , (fePi), with (di)- 
symmetric Pi for i = 1, . . . , n, and where C is a or (di (^-symmetric. A computation similar 
to that of eqn. (12. 4p yields the assertion. 

6. Staircase polygons are also characterised by the property that they contain the lower left and 
the upper right corner of their smallest bounding rectangles. So the only staircase polygons with 
(h)- or (v)- symmetry are rectangles. / maps a rectangle Q either to a single rectangle, or to 
sequences of vertical (horizontal) edges, if the width (height) of Q is 1. This results in the above 
equation. 

7. The only admissible polygons are squares. For a given half-perimeter, there is exactly one 
square. If n > 1, the function / maps a square of half-perimeter 2n to the square of half-perimeter 
2n — 2, and it maps the unit square to pt. We obtain the claimed equation. 

□ 

Remark 2.2. Equations of the above form appear in different contexts. Examples are classes 
of directed lattice paths, counted by length and area under the path [32], or classes of simply 
generated trees |31) . counted by number of vertices and internal path length. This is due to 
combinatorial bijections between these classes, which we will partly review in Section [5] In the 



5 



context of polygon models, equations appear for Class 1 in [3] and for Class 6 in [37], while Class 
7 is trivial. Solutions of some equations may be given in explicit form, compare J5[ 134] . 



3 Area limit laws 

In this section, we derive the limiting area laws for the various symmetry subclasses, in the 
uniform fixed perimeter ensemble. This will be achieved by an application of the moment method 
[3j Sec. 30]. Such an approach has been used previously [42] [43j [32l [12] in similar contexts, 
using some involved computations. We will follow a streamlined version, based on the method 
of dominant balance [38] . which finally allows to obtain the limit distribution by a mechanical 
calculation. In order to give a self-contained description of the method, we will treat the two cases 
V and S in detail, and then indicate the analogous arguments for the remaining subclasses. 



3.1 Limit law for V 

A g-difference equation for the half-perimeter and area generating function P(x,q) of all staircase 
polygons was derived in Proposition 12.11 For q = 1 , the resulting quadratic equation describes 
the generating function P{x, 1) of staircase polygons, counted by half-perimeter. The relevant 
solution is 

P(x,l) = i-ivT=4£+i(l-4aO. (3.1) 

We are interested in the distribution of area within a uniform ensemble where, for fixed perime- 
ter 2m, each polygon has the same probability of occurrence. We introduce the discrete random 
variable X m of area by 

where [u k ]f(u) denotes the coefficient of order k in the formal power series f(u). In the following, 
we will asymptotically analyse the moments of X m . The answer can be expressed in terms of the 
Airy distribution, see |15[ I24( 125] for a discussion of its properties. 

Definition 3.1. A random variable Y is Airy distributed 17 5^ if 

E[y fc ] = r( 7o ) (f> k 

kl r(7fe)0 o ' 

where 7^ = 3fc/2 — 1/2, and where T(z) is the Gamma function. The numbers 4>k satisfy for k £ N 
the quadratic recursion 

1 k 

2=0 

with initial condition 4>o = — 1- 

Remark 3.2. i) In the sequel, we shall make frequent use of Carleman's condition: A sequence 
of moments {M m } mS N with the property ~^2 k {M2k)~ 1 ^ 2k ^ = 00 defines a unique random variable 
X with moments M m , cf. [15] . 

ii) This implies in particular, that Y is uniquely determined by the above moment sequence. 
Explicit expressions can be given for its moments, its moment generating function, and its density. 
The name relates to the asymptotic expansion 

A logAioo ~ E(- 1 ) fc ^ s " 7fe ( s ^ °°). 

fe>0 

where Ai(x) = \ / °°cos(t 3 /3 + tx)dt is the Airy function. The Airy distribution appears in a 
variety of contexts. In particular, the random variable Y/y/8 describes the Brownian excursion 
area. 



G 



We can now state the following result. 

Theorem 3.3. For staircase polygons of half-perimeter m, the area random variables X m eqn. (|3.2| . 
appropriately normalised, converge in distribution, 

X m d Y 

where Y is Airy distributed. We also have moment convergence. 

Remark 3.4. The previous theorem is a special case of [12l Thm. 3.1] and [381 Thm. 1.5]. In 
[12j . a limit distribution result is stated for certain algebraic g-difference equations, together with 
arguments of a proof using the moment method. In [38j . a general multivariate limit distribution 
result is proved for certain q-functional equations, using the moment method and the method of 
dominant balance. The corresponding argument for staircase polygons, using the moment method 
and the method of dominant balance, is sketched in [55] , 

For pedagogic purposes, we will give a proof of the above result using the moment method and 
the method of dominant balance. In particular, we show that the moments of X m , appropriately 
normalised, converge to those of Y. Since Y is uniquely determined by its moments by Carle- 
man's condition, this implies that the sequence {X m } m£ jq converges, after normalisation, to Y in 
distribution and for moments, compare JS[ Thm. 4.5.5]. 

We will analyse the asymptotic behaviour of the moments in terms of the singular behaviour 
of the associated factorial moment generating functions which are, for (fc, () 6 Nq, defined by 



Qk+l 



9=1 



In particular, Po t o(x) — P(x, 1) is the generating function of staircase polygons, counted by half- 
perimeter. Note that these quantities exist as formal power series and have the same radius of 
convergence as P(x, 1). This is due to the fact that the area of a polygon is bounded by the 
square of its perimeter, resulting in [x m ] P(x, q) being a polynomial in q. The name results from 
the identity 

( A mjfc - r m1 5 — r^r, (3.3J 
[x m \Pou{x) 

where (a)k = a ■ (a — 1) • . . . • (a — k + 1) is the lower factorial. The factorial moment generating 
functions Pk,i{x) turn out to be algebraic. Explicit expressions may be obtained recursively from 
the functional equation eqn. (|2.ip , by implicit differentiation w.r.t. x and q. 

We will study the singular behaviour of the factorial moment generating functions from their 
defining functional equation, and will then infer the asymptotic behaviour of the moments. This 
will be achieved in three steps. We will first prove the existence of a certain local expansion for 
each factorial moment generating function about its singularity, by an application of the chain 
rule (or Faa di Bruno's formula). Then we will provide an explicit expression for the leading 
term in the expansion, by an application of the method of dominant balance. This will finally be 
analysed in order to obtain the asymptotic behaviour of the corresponding moment, by methods 
from singularity analysis of generating functions. 

We note that steps one and two are usually performed simultaneously, the corresponding 
method being nicknamed moment pumping [T7]. Our two-step approach uses an exponent guess, 
which might be obtained from an analysis of the first few factorial moment generating functions. 
It is then shown that the guessed exponent value is an upper bound on the true exponent, by 
an application of Faa di Bruno's formula. The corresponding calculation is simpler to perform 
than the usual asymptotic analysis of the functional equation. Once an exponent bound has 
been established, the method of dominant balance can be applied. It yields a recursion for the 
coefficients of the leading singular term in the factorial moment generating functions. If the 
recursion reveals non-zero coefficients, this proves that the exponent bound is actually an equality, 
thereby verifying the initial guess. 
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The first step of our method is summarised by the following lemma. For its statement, recall 
that a function f(u) is A-regular [14] if it is analytic in the indented disc A = A(tt c ) = {u : \u\ < 
u c + rj, |arg(u — u c )\ > (/)} for some real numbers u c > 0, n > and <fi, where < <f> < 7r /2- Note 
that u c ^ A, where we employ the convention arg(0) = 0. The set of A-regular functions is closed 
under addition, multiplication, differentiation, and integration. Moreover, if f(u) ^ in A, then 
1/ f(u) exists in A and is A-regular. 

Lemma 3.5. For (k, I) S Nq, the power series Pk,i(x) has radius of convergence 1/4 and is A(l/4)- 
regular. It has a locally convergent expansion about x = 1/4, as in eqn. (|3.1|) for (k,l) = (0,0), 
and for (k, I) ^ (0, 0) of the form 

oo ^ 

Pk.l (X) = (1 _ 4a . )3fc /2+ i -r/2-l/2 ■ ( 3 - 4 ) 

r— ^ ^ 

Remark 3.6. i) The exponent 3k/2 in eqn. (|3.4j) might be guessed from the asymptotic behaviour 
of the mean area E m [X m ] ~ Arr?! 2 of a random polygon. The mean area is obtained from 
Pa,o{x) and P±fi{x), which might be easily extracted from the g-difference equation. Note that 
the coefficients dk,i,o might attain zero values at this stage. The recursion eqn. (|3.10p given below 
however implies that all of them are non-zero. 

ii) The reasoning in the following proof may be used to show that all series Pk,i{%) are algebraic. 
Hi) For (r 2 ) -symmetric polygons, our proof below will use properties of the derivatives 

These functions have all radius of convergence 1/2, are A(l/2)-regular, and have the same type of 
expansion as the functions Pk,i{x) of the previous lemma. This may be inferred from the previous 
lemma by the chain rule or, more formally, by an application of Faa di Bruno's formula [8j. 

Proof. By the argument given above, it is seen that all functions Pk,i{x) have the same radius 
of convergence. The statement of the theorem is true for Po,o(aO = P(x,l), as follows from the 
explicit expression eqn. (13. ip . For the general case, we argue by induction on (k, I), using the total 
order < defined by 

(r, s) < (k, I) r + s < k + I V (r + s = k + 1 A r < k) , 
chosen to be compatible with the combinatorics of derivatives. Define 

H(x, q) := P(x, q)(l - 2xq - P(xq, q)) - x 2 q, 
compare eqn. (|2.ip . Fix (fc, I) t> (0, 0). An application of Leibniz' rule yields 
d k+l , . ,s fk\fl\ d k+l - r ~ 6 



{H^q) +X \)= \ [)[)——— A P {x , q)) 



dq k dx l ' \r J \s J dq k r dx l 



(0,0)<(r,s)<(fc,0 



(3.5) 

Qr+s 

(l-2xq-P(xq,q)) 



dq r dx s 

In fact, terms corresponding to indices (r, s) < (fc, I) with s > I or r > k are zero. For the second 
derivative on the r.h.s. of eqn. (|3.5p . note that by the chain rule 

^p(.«»=t(;)' r - i (s^ p )< i «i' 

Taking further derivatives w.r.t. x, we may write 

(P( X q,q)) = q^ -—-p)(q Xl q) 



dq r dx s ' \dq r dx s 

+ L {wd^ p {xq > q) - WiAx ' q) > 
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for polynomials Wi t j[x, q) in x and g, which satisfy Wij{x, q) = if i < r. By inserting eqn. (|3.6p 
into eqn. (|3 . 5|) and setting g = 1, one observes that only the (0,0) and the (k,l) summand in 
eqn. Q3.5P contribute terms with Pkj(x). The terms involving Pk,i{x) sum up to 



P k ,i(x) (1 - 2z - 2P ,o(z)) = v 7 ! - 4^P M (a;), 

where we used eqn. (|3.1J . Now the claimed A(l/4)-regularity of Pk,i{x) follows from the induction 
hypothesis, by the closure properties of A-regular functions. For the particular singular expansion 
eqn. (|3.4[) note that, by induction hypothesis, each of the remaining terms in the summation in 
eqn. (|3.5p has an expansion eqn. (|3.4[) . Hence, the most singular exponent is bounded by 

|(fc-r) + (i--)-|) + (|r + --|) =§* + /-!. 

We conclude that the leading singular exponent of Pk,i(x) is at most 3k/2 + I — 1/2, which yields 
the desired bound, and thus the remaining assertion of the theorem. 

□ 



The second and third step of our method yield a proof of Theorem 13.3 



Proof of Theorem \3.3l We apply the method of dominant balance [38] in order to obtain the limit 
distribution of area. Its idea consists in first replacing the factorial moment generating functions, 
which appear in the formal expansion of P(x,q) about q = 1, by their singular expansion of 
Lemma l3.5[ and then in studying the equation implied by the g-difference equation eqn. (|2.1[) . We 
may thus write 

P(X, q) = \ + (1 - q)^ F ( -L^, (1 - g)^ , (3.7) 

where F(s,e) — ^2 F r (s)e r is a formal power series in e with coefficients F r (s) being formal 
Laurent series in s 1 ' 2 . The series 

k 

is some generating function for the leading coefficients of Pf~o(x). The coefficients fk '■= rffc,o,o/^!> 
in turn, determine the asymptotic form of the factorial moments E[(X TO )fc] in eqn. (|3.3p . as we 
will see below. We will use the g-difference equation eqn. (|2.ip to derive a defining equation for 
Fq(s). This will lead to a simple quadratic recursion for the numbers fk- Use the above form of 
P(x, q) in the g-difference equation, introduce Ax = 1 — se 2 , q = 1 — e 3 , and expand the functional 
equation to second order in e. This yields a Riccati equation for the generating function Fq(s), 

^-F (s) + AF (s) 2 -a = 0. (3.9) 
as 

On the level of coefficients of F (s), we obtain the recursion 

fe 

7wA-i + 4^/i/w = 0, (3.10) 
;=o 

with initial condition fo = do,o,o = ~ 1/2. We infer from the definition of the Airy distribution 
that fk = 2~ 2k ~ 1 (f>k. In particular, all coefficients fk are non-zero. Noting that the functions 
Pk,o( x ) are A(l/4)-regular, we thus get by an application of the transfer lemma [THl Thm. 1] for 
the factorial moments of X m the asymptotic form 

EmlWk] _ 1 [x m ]Pk,o(*) 1 [^ m ]4,o,o(l - 4z)-( 3fe / 2 -V2) 



fc! fc! [x m ]P 0fi {x) k\ [x m ](io,o,o(l ~ 4a;) 1 / 



2 



fo r(3fc/2 - 1/2) 0o r( 7fc ) V 4 ; 1 
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The previous estimate also shows that the factorial moment E m [(X m )fe] is asymptotically equal 
to the ordinary moment E TO [(X TO ) fc ]. It follows with Thm. 4.5.5] that the sequence of random 
variables {4m -3 / 2 X m } mgN converges in distribution to Y, where Y is Airy distributed. The above 
reasoning also implies moment convergence. 

□ 



Remark 3.7. If we expand the functional equation to higher order in the above example, we 
obtain at order e r+2 a linear differential equation for the function F r (s), which is the generating 
function for the numbers dk,o,r in the expansion eqn. (|3.4p . compare [37j . So we can mechanically 
obtain corrections to the asymptotic behaviour of the factorial moment generating functions, and 
hence to the moments of the limit distribution. 



3.2 Limit law for S 

The above strategy can also be followed in order to study the area law for the class of (r 2 )- 
symmetric staircase polygons. The result can be expressed in terms of the distribution of area of 
the Brownian meander, see [43l Thms. 2,3] and the review [24] . 

Definition 3.8. The random variable Z of area of the Brownian meander is given by 

E[Z k ] = T(a ) uj k 1 
k\ T(a k ) ojq 2 fc / 2 ' 
where a k = 3k/ 2 + 1/2. The numbers u>k satisfy for k € N the quadratic recursion 

k 

a k -iuJk-i +^0;2~W-i = 0, (3-11) 

with initial condition loq — 1, where the numbers <pk appear in the Airy distribution. 

Remark 3.9. By Carleman's condition, the random variable Z is uniquely determined by its mo- 
ments, and explicit expressions are known for the moment generating function and the distribution 
function. 

We are particularly interested in the derivatives 

where (k, I) £ N 2 ,. As above, these series exist as formal power series and have the same radius of 
convergence. We have the following lemma. 

Lemma 3.10. For (k,l) £ N 2 ,, the power series Sk,i(x) has radius of convergence 1/2 and is 
A(l/2) -regular. It has a locally convergent expansion about x = 1/2 of the form 

SfcJ.r 



Sk,i(x) ^2 . _ 2 ., 2+ , 2l , 2 • 

r>0 v ' 

Remark 3.11. The following proof can be used to show that all series Sk,i{x) are algebraic. 



Proof. The proof is analogous to that of Lemma 13.51 Elementary estimates show that all series 
Sk,i(x) have the same radius of convergence. Setting q = 1 in eqn. (|2.2p . solving for So t o(x) and 
expanding about x — 1/2 yields the assertion for (k,l) = (0,0). We argue by induction on (k,l), 
using the total order <]. Fix (k,l) > (0,0). Differentiating eqn. (|2.2[) with Leibniz' Rule gives 

d k+l „, , fk\fl\ d r+s fP(x 2 ,q 2 ) 

r S{x,q)- x ' 



k dx l ' ^— ' \r J \s I dq r dx s V x 2 q 

(o,o)<(r,»)<(fc,i) v 7 v 7 v - ( 3 _ 12 ) 

Qk+l—r—s 

(1 + 2xq + S(xq,q)) . 



Qgk-rQ x l- 
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We argue as in the proof of Lemma [3. 51 that only the (0,0) summand on the right hand side of 
eqn. (|2.2p contributes (fc, I) derivatives of S, and that all other derivatives of S of order (r, s) 
satisfy (r,s) < (k,l). Setting q = 1 in (|2.2[) and collecting all terms involving Sk,i(x) on the left 
hand side gives 



p o,o( x ) \ a ( \ V^) n , n , 
—2—\S ktl {x) = —2—(l + 2x)+ ^ 

/ ' (*•>«)< O,0 
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r+s 



dq r dx s 



P(x 2 ,q 2 
x 2 q 



9=1 



(3.13) 



hr t s{x) ~\~ ^ ' O-i.jSij(x) 
(i,j)<(k,l) 



where the h r ^ s {x) are (at most linear) polynomials, and the a^j are some real coefficients. Note 
also that the terms 



0rH 



dq r dx s 



P(x 2 ,q 2 ) 
x 2 q 



9=1 



E 



Pi A*) 



x 2+r-iql + s—j 



are A(l/2)-regular, with an expansion about x = 1/2 having the same exponents as in eqn. 
see the remark following Lemma [33] We thus get A(l/2)-regularity of Sk.i(x) by induction, and 
by the closure properties of A-regular functions. For the particular expansion, note that the right 
hand side has a locally convergent expansion about 1/2 with most singular exponent 3A:/2 + I, as 

the factor 



dx* ( 



9=1 



has an expansion with most singular exponent 3r/2 + s— 1/2, and 

the inner sum has by induction an expansion with most singular exponent at most 3(fc — r)/2 + 
(I — s) + 1/2. The first factor on the left hand side has a locally convergent expansion about 1/2 
starting with 

PoA*) X 



1 - 



= -2v / 2v / l - 2x + 0(1 - 2a 



1/2). 



Solving for Sk,i(x) yields the desired expansion. 



Theorem 3.12. The random variables X 



(syr. 



□ 



of area of (r 2 ) -symmetric staircase polygons of 



half-perimeter ra, appropriately normalised, converge in distribution, 



v{sym) ry 



m 



3/2 



— (m 

2 V 



oo), 



where Z is the meander area random variable. We also have moment convergence. 

Proof. We apply the method of dominant balance. According to Lemma I3.10[ the generating 
function S(x, q) may be expressed as 



S(x,q) 



G 



1 - 2x 



',(1-9) 



1/3 



(1-?)V3 \(l-q)V3 

where G(s, e) = Yl G r (s)e r is a formal power series in e and s -1 / 2 , and 

Sfc,o,o (-!) fc 



(3.14) 



G(s,0) = G (s)=J2- 



fcl s 3fe/2+l/2 



is a generating function for the leading coefficients in the singular expansions of the functions 
Sk,o(x). The functional equation eqn. (|2 . 2[) induces a recursion on the numbers gk '■= Sk,ofi/kl, 
which determines the limit distribution, as we will see below. We insert eqn. (|3.14j) together with 
eqn. ()3.7j) into the functional equation, introduce q = 1 — e 3 and 2x — 1 — se 2 , and expand the 
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functional equation to order zero in e. This gives the linear inhomogeneous first order differential 
equation 

^Go(s) + 4 ■ 2 1 / 3 F Q (2 1 / 3 s) G (a) + 2 = 0, (3.15) 
where Fq(s) is given by eqn. (|3.8p . On the level of coefficients, we have the recursion 

k 

ak-x9k-i + 2 ~ l/2+5/2 fi9k-i = 0, (3.16) 

1=0 

with the initial condition go — so,o,o = 2 -1 / 2 . If we set 



9k 



23fc/2+l/2 ' 



then the above recursion is identical to that occurring in the definition of the meander distribution. 
In particular, all numbers gk are non-zero. Since the functions Sk,o(x) are A(l/2)-regular, we may 
use the transfer lemma [TH1 Thm. 1] to infer for the moments of li 59 " 1 ' the asymptotic form 

E m [(X^ m) ) fc ] _ 1 [x m ]S k , (x) 1 [s m ]s M , (l - 2x)-( 3fc / 2 +V 2 ) 
k\ k\ [x m ]S ,o{x) ~ k\ [x m ]s , ,o(l - 2a;)- 1 /2 

5o T(3fc/2 + 1/2) 2 fe w r(a fc )2 fc / 2 V ; 

The last term is, up to the factor m 3fc//2 , the &;-th moment of Z/2, where Z is the meander area 
variable. The previous estimate shows that the factorial moments ¥. m [{Xm Vrr ^)k] are asymptot- 
ically equal to the ordinary moments K m [(Xm V ) k ]- It follows with [SJ Thm. 4.5.5] that the 
sequence of random variables j 2m" 3 / 2 lm ! '™'' 1 converges in distribution to Z, where Z is 
distributed as the meander area. We also have moment convergence. 

□ 

Remark 3.13. As for the full class of staircase polygons, corrections to the asymptotic behaviour 
of the factorial moment generating functions can be mechanically obtained also for this example, 
by expanding the corresonding functional equation to higher orders in e. 

3.3 Limit law for (di)-symmetric polygons 

These polygons always have even half-perimeter. In order to derive an area limit law, we thus 
restrict to area random variables of quarter-perimeter. Note that we have D\{x, q) — _Di(x 1 / 2 , q) 
for the generating function of the class (di)-symmetric polygons, counted by quarter-perimeter 
and area. The functional equation for D\(x,q) induces a similar one for D\(x,q). Their factorial 
moment generating functions all have radius of convergence 1/4, and a statement as in Lemma 
13.51 can be formulated and proven almost verbatim for D\(x, q). The method of dominant balance 
then yields a generating function for the leading coefficients in the singular expansions, a defining 
equation similar to eqn. (|3.9p . and a recursion similar to eqn. (|3. 10[) . We have the following result. 

Theorem 3.14. The area random variables X m of (di)- symmetric staircase polygons, indexed by 
quarter-perimeter m and scaled by m~ 3 / 2 , converge in distribution to a random variable Y, which 
is Airy distributed. We also have moment convergence. 

3.4 Limit law for (c^-symmetric polygons 

In [26] . a combinatorial bijection between (c^-symmetric polygons and (r 2 )-symmetric polygons 
with even half-perimeter is described: cut a ((^-symmetric polygon along the line of reflection, 
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flip its upper right part, and glue the two parts together along the cut. So Theorem l3 . 1 21 translates 
to the (c?2)-case. 

Alternatively, one may apply the methods of Section T3. 21 together with modifications similar to 
those of Section \3. 31 to the quarter-perimeter and area generating function D2{x, q) = D2(x 1 ^ 2 , q). 
Lemma [3.101 holds in this case, with 1/2 replaced by 1/4, and the method of dominant balance 
yields results similar to eqn. (|3.15p and eqn. (|3. 16|) . 

Theorem 3.15. The area random variables X m of (d^)- symmetric staircase polygons, indexed by 
quarter-perimeter m and scaled by (2m)~ 3 / 2 , converge in distribution to a random variable Z/2, 
where Z is the meander area random variable. We also have moment convergence. 

3.5 Limit law for (d\, ^-symmetric polygons 

In this symmetry class, every polygon has even half-perimeter. So we define Dviix, q) = D^x 1 ^ 2 , q) 
as above, and obtain from the functional equation for D\2(x,q) one for D\2{x,q), which involves 
Di(x,q), resembling eqn. 1)2.2|) . 

It can be argued, as in the proof Lemma 13.101 that all factorial moment generating functions 

„k — 

4-rrDi2{x,q) have radius of convergence 1/2, with singularities at ±1/2, where the leading 
q ^T 1 

singular behaviour of the coefficients is determined by the singularity at 1/2. We can apply the 
methods of Section [3~2l with the modifications of Section [3~3l This yields the following result. 

Theorem 3.16. The sequence of area random variables X m of (d\, d2) -symmetric staircase poly- 
gons, indexed by quarter-perimeter m and scaled by m~ 3 / 2 , converges in distribution to 2Z, where 
Z is the meander area variable. We also have moment convergence. 

3.6 Limit law for (r)-symmetric polygons 

The class of staircase polygons with (r)-symmetry is the class of squares. These may be counted 
by quarter-perimeter m. Since for given m there is exactly one square, they have, after scaling 
by m~ 2 , a concentrated limit distribution 5(x — 1). This result can also be obtained from the 
q-difference equation in Proposition l2.il 

3.7 Limit law for (h,v)- ((h)-, (v)-) symmetric polygons 

The class of staircase polygons with (h, w)-symmetry (or with (h)- or (w)-symmetry) is the class 
of rectangles. These have been discussed in |40j . The fc-th moments of the area random variable 
X m , with m half-perimeter, cf. eqn. (13. 2p . are given explicitly by 

^ 1 f 1 (kU 2 

E[X* ] = J2 W m - *)) fc ^3T ~ ^ J a " x ^ dx = (2fcTT)! m2fe {m ~+ 00)1 

where we used a Riemann sum approximation. Consider the normalised random variable 

X m = 4X m /m 2 . 

The moments of X m converge as m — > 00, and the limit sequence Mj~ — limm^oo E m [X^] sat- 
isfies Carleman's condition and hence defines a unique random variable with moments Mk- The 
corresponding distribution is the beta distribution /3i.i/2- We arrive at the following result. 

Theorem 3.17. The sequence X m — AX m /m 2 of area random variables of rectangles, with half- 
perimeter m scaled by 4/m 2 , converges in distribution to a 0x t x/ '2- distributed random variable. We 
also have moment convergence. 

One may also obtain this result by manipulating the associated (/-difference equation, sec [40 . 
Expansions of the factorial moment generating functions about their singularity at x = 1 can be 
derived, and bounds for their most singular exponent can be given. The method of dominant 
balance can then be applied to obtain the leading singular coefficient of these expansions. 
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4 Limit law for orbit counts 



Let TL be a subgroup of V4. By the Lemma of Burnside, the half-perimeter and area generating 
function Pu{x,q) of orbit counts w.r.t. TL is given by 26J 

Pn(x,q) = p FiK(g)(x,q) = (P(x,q)+R(x,q)) , 

where |?i | denotes the cardinality of TL, and where Fix(g) C V is the subclass of staircase polygons, 
which are fixed under g £ 2?4, with half-perimeter and area generating function Pp; x ( s )(x, q). The 
series P(x, q) is the full staircase polygon half-perimeter and area generating function, and R(x, q) 
is the sum of generating functions of staircase polygons which are fixed under g £ TL, where j^c. 
Let P(x, 1) = J2mP mX " 1 an d R{x, 1) = S m r m a;m - Due to the previous discussion, see also [26l 
Prop. 14], the number of polygons fixed by some non-trivial symmetry grows subexponentially 
w.r.t. the total number of polygons. This implies 

-----^0 (m->oo), (4.1) 

Prn 

for any real number a. As a consequence, area limit distributions for orbit counts coincide with 
those for the full class of staircase polygons. 

Theorem 4.1. Let TL be a subgroup of T>±. Then the area limit law of the class V /TL coincides 
with that of V . 

Proof. We show that both classes have asymptotically the same area moments. Since in all 
examples the limit distribution is uniquely determined by its moments, the claim follows. 

Note that, for polygons of half-perimeter m, their area n satisfies 1 < n < m 2 . Let P(x 1 q) — 
^2,Pm,nX m q n and R(x, q) =J2 r m, n x m q n . By eqn. (|4. 1|) . we have for fceNo 



rn r 



< > (m — > 00). 



This implies for the coefficients of the moment generating functions the asymptotic estimate 

\ Xm ] ( IJT ) P n(x, q)\ q=1 = 7^77 nk (jPrn,n + r m ,„) 



1 / rl \ k 

(m — > 00). 



^E"W = pM 



We conclude that both classes have asymptotically the same area moments. □ 



5 Staircase polygons, Dyck paths, and Brownian excursions 

We briefly explain how some of our results could be alternatively obtained from the bijections 
described in |26j . which set symmetry classes of staircase polygons in one-to-one correspondence 
to symmetry classes of Dyck paths. Random Dyck paths, in turn, are related to corresponding 
stochastic objects such as the Brownian excursion and the Brownian meander. This allows to infer 
area limit distributions for some polygon symmetry classes from distributions of certain Brownian 
excursion and meander functionals. 

A combinatorial bijection between staircase polygons of perimeter 2m + 2 and Dyck paths 
of length 2m has been described by Delest and Viennot [5]. Within that bijection, the area 
of a staircase polygon corresponds to the sum of the peak heights of a Dyck path. In a uniform 
ensemble, the sequence of random Dyck paths w.r.t. half-length yields, after suitable normalisation, 
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a sequence of stochastic processes on (C[0, 1], || • ||oo) with the Borel c-algebra, which converges 
in distribution to the standard Brownian excursion [HH]. Also, for certain continuous functionals 
on suitably normalised Dyck paths, including the area and the sum of peak heights, we have 
convergence in distribution and moment convergence to the Brownian excursion area, see [55, 20J 
and pT[ Thm. 9]. Using stochastic techniques, the distribution of the excursion area has been 
initially analysed in STJ [29] , compare the historical remarks in [33] . Derivations by discrete 
methods are summarised in |15j . 

The above bijection, restricted to r 2 -symmetric polygons, yields symmetric Dyck paths, which 
decompose in two identical discrete meanders [26] . Again, it is known that the sequence of random 
discrete meanders w.r.t. half-length, appropriately normalised, converges to the Brownian meander 
23J . Together with the convergence theorem for continuous functionals of polynomial growth [111 
Thm. 9], one may conclude that also discrete meander functionals, such as area and sum of peak 
heights, converge to the Brownian meander area. The distribution of the Brownian meander area 
was derived in [53] from a ^-difference equation similar to eqn. (|2.2p . It can also be obtained 
from known results for the corresponding Brownian motion and bridges distribution, since there 
is a relation between the double Laplace transforms of the three distributions [35] [24]. The 
same considerations hold for G^-symmetric polygons, which are in one-to-one correspondence to 
r 2 -symmetric polygons [26] . 

Staircase polygons with c?i-symmetry are in bijection with pairs of identical Dyck paths, which 
is seen by cutting a polygon along its positive diagonal [26j . Here, the polygon area corresponds 
to twice the Dyck path area, hence the area limit distribution is given by that of the Brownian 
excursion area, by the convergence result mentioned above. A similar result holds for (di,^)- 
symmetric polygons. Here, a combinatorial bijection between polygons and discrete meanders is 
known [26] , where the polygon area corresponds to four times the meander area. 

The classes of rectangles and squares lead to Dyck paths with an initial sequence of up steps, 
followed by an alternating sequence of peaks and valleys, ending in a terminal sequence of down 
steps. As seen above, these classes are treated by elementary methods. 

6 Conclusions 

We analysed the symmetry subclasses of staircase polygons on the square lattice. Exploiting a 
simple decomposition for staircase polygons [39] , we obtained the area limit laws in the uniform 
fixed perimeter ensembles. This extends and completes previous results [26] , As expected, orbit 
counts with respect to different symmetry subgroups always lead to an Airy distribution. The 
enumeration of polygons fixed under a given symmetry group leads to a variety of area limit 
distributions, such as a concentrated distribution, the ^ ^-distribution, the Airy-distribution, 
or the Brownian meander area distribution. 

As described in Section [5] the latter two results can also be obtained from the connection 
to Brownian motion and the Brownian meander. Our independent discrete approach uses only 
elementary methods from probability and singularity analysis of generating functions. Moreover, 
it may be used in order to analyse corrections to the limiting behaviour, which cannot easily be 
obtained by stochastic methods. 

One may also study the analogous problem of perimeter limit laws in a uniform ensemble 
where, for fixed area, every polygon occurs with the same probability. For the class of staircase 
polygons, the associated centred and normalised random variable is asymptotically Gaussian [18[ 
Prop. 9.11], and the same result is expected to hold for the symmetry subclasses, apart from 
squares. Also, limit laws in non-uniform ensembles and for other counting parameters may be 
studied, compare [40] . 

The above methods may be applied to extract limit laws related to symmetry subclasses of 
other polygon classes. In particular, the classes of convex polygons on the square [27] and on the 
hexagonal [5T] lattices may be studied. 

With respect to symmetry subclasses of self-avoiding polygons, exact enumeration studies may 
be carried out. It would be interesting to numerically analyse moments, in order to test conjectures 
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for limit distributions, compare [361 140) . 
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